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ABSTRACT 

Objectives: Evidence on the association between the 
adverse socioeconomic characteristics of residential area 
and mortality is mixed. We examined whether the choice 
of spatial unit is critical in detecting this association. 
Design: Register-linkage study. 
Setting: Data were from the Finnish Public Sector study's 
register cohort. 

Participants: The place of residence of 146 600 cohort 
participants was linked to map grids and administrative 
areas, and they were followed up for mortality from 2000 
to 2011. Residential area socioeconomic deprivation and 
household crowding were aggregated into five alternative 
areas based on map grids (250x250 m, 1x1 km and 
10x10 km squares), and administrative borders (zip-code 
area and town). 

Primary and secondary outcome measures: 

All-cause mortality. 

Results: For the 250x250 m area, mortality risk 
increased with increasing socioeconomic deprivation (HR 
for top vs bottom quintile 1.36, 95% C1 1.21 to 1.52). 
This association was either weaker or missing when 
broader spatial units were used. For household crowding, 
excess mortality was observed across all spatial units, the 
HRs ranging from 1.14 (95% C1 1.03 to 1.25) for zip 
code, and 1.21 (95% C1 1.11 to 1.31) for 250x250 m 
areas to 1 .28 (95% C1 1 .1 0 to 1 .50) for 1 0x1 0 km 
areas. 

Conclusions: Variation in spatial units for analysis is a 
source of heterogeneity in observed associations between 
residential area characteristics and risk of death. 



ARTICLE SUMMARY 



Article focus 

■ There is no strong consensus on which spatial 
units are best for determining the health effects 
of residential areas. 

■ Few studies have been able to compare area-level 
socioeconomic effects using several alternative 
spatial units. 

Key messages 

■ Data on residential area socioeconomic depriv- 
ation and household crowding were aggregated 
into five alternative areas based on map grids 
(250x250 m, 1x1 km and 10x10 km squares), 
and administrative borders (zip-code area and 
town/city). 

■ High areal socioeconomic deprivation and house- 
hold crowding, as aggregated into the smallest of 
the five spatial units, 250x250 m square, were 
associated with increased mortality. For household 
crowding, excess mortality risk was also observed 
using the other spatial units. 

■ These data show that aggregating data in differ- 
ent ways leads to different results in the analyses 
of the associations between residential area char- 
acteristics and risk of death. 

Strengths and limitations of this study 

■ Individual socioeconomic variables were adequately 
controlled for. 

■ As the study population consisted of Finnish public 
sector employees, the generalisability of the results 
needs to be confirmed in other studies. 
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INTRODUCTION 

Evidence that the adverse socioeconomic char- 
acteristics of residential areas are risk factors 
for all-cause mortality is mixed, comprising 
both positive 1-21 and null findings. 6 14 22 In 
these studies the spatial unit to which area data 
has been aggregated has varied considerably 



and is a possible source of inconsistencies, a 
feature known as the Modifiable Area Unit 



Problem (MAUP). 
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Some investigations 



have aggregated area characteristics to the level 
of states 25 towns 14 22 zip-code areas 11 21 2b 27 
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units.' ' Towns and other large administrative 

units can capture differences in the provision of commu- 
nity health and welfare services, but smaller spatial units, 
such as zip codes, may cover local variability in people's 
social environments as well as 'local health-related 
cultures' that may also contribute to mortality differences 
between areas. 

Prior research comparing health effects by spatial units 
has suggested that no differences exist between spatial mea- 
sures 11 16 27 29 or that the smaller ones provide stronger 

4 13 18 28 

effect estimates. ' However, few studies have systemat- 
ically examined this issue across different area characteris- 
tics and various spatial units within a single analytic setting 
and adequately adjusting for individual socioeconomic vari- 
ables. We sought to undertake such a study by comparing 
five different spatial units (towns, zip-code areas and 
map-grid squares of 250x250 m, lxl km and 10x10 km) in 
relation to two widely used socioeconomic area characteris- 
tics, deprivation and household crowding. 



METHODS 

Study design and population 

The Finnish Public Sector study cohort consists of employ- 
ees working for ten municipalities and six hospital districts 
in Finland. All men and women employed in these organi- 
sations for more than 6 months in any year between 1991 
and 2005, and from the full spectrum of socioeconomic 
groups were eligible (n=151 901). Owing to the nature of 
public sector jobs in Finland (nurses, teachers, etc) most of 
the study participants were women. For this study, we 
selected those cohort members who were alive and aged 
18-65 years at the beginning of the follow-up, which was 
the date on which the participant began his/her first 
employment contract in the target organisations between 1 
January 2000 and 1 January 2005 (for those contracted 
before 2000, the start date was 1 January 2000) . Thus, of 
the included participants 95.5% were employed at the 
beginning of the follow-up. Participants were followed until 
the end of December 2011, a move abroad, or death, 
whichever came first. The global positioning system 
(GPS)-coordinates of the residential buildings of 146 831 
participants were obtained from the population informa- 
tion system of the Population Register Center using per- 
sonal identification codes. The centre's data on nearly 
three million residences is maintained and checked in 
close cooperation with municipal building supervision 
authorities and local register offices. 31 The outcome was all- 
cause mortality, including deaths from diseases and exter- 
nal causes. The dates of death were obtained from Statistics 
Finland on the basis of personal identification codes. The 
Ethics Committee of the Hospital District of Helsinki and 
Uusimaa approved the study. 



Area characteristics and spatial units 

The spatial units used were map-grid-based squares 
of 250x250 m, lxl km and 10x10 km and administrative 
measures of town and zip-code areas. 32 The smallest units 



represent areas in which people communicate with their 
neighbours, and conduct their businesses by foot. The 
zip-code area is a unit used in prior literature 4 11 21 and is 
based on a defined postal area, usually larger than 1 km 
square, but smaller than a town. Towns are thus assumed 
to form boundaries within which people conduct most of 
their daily activities, and larger grids, in the 10x10 km 
scale, represent units possibly comparable to administra- 
tive units. All these units cover the whole Finland. 

The participants were linked to the map squares by using 
the GPS-coordinates of their residential buildings and to 
their administrative areas by their postal addresses. For 
each spatial unit, information on deprivation and house- 
hold crowding was calculated by Statistics Finland on the 
basis of data from the population register, registers of the 
Finnish Tax Administration, and Statistic Finland's employ- 
ment register. In all these registers the total population res- 
iding in Finland at the end of data collection year served as 
the universe. For each spatial unit, we defined an index of 
socioeconomic deprivation using information on median 
income (median household income in the area logarith- 
mically transformed and then coded as additive inverse in 
order to obtain higher values for greater deprivation) , edu- 
cation attainment (proportion of those aged >18 whose 
highest education level was elementary school), and 
unemployment rate (unemployed people belonging to the 
labour force/ total labour force). These are standard vari- 
ables, used either separately or joindy, to characterise the 
disadvantage and deprivation. 33 For each of the three indi- 
cators, we derived a standardised z-score (mean=0, SD=1). 
The index of socioeconomic deprivation was then calcu- 
lated by taking the mean value across all z-scores 34 when 
the z-score for at least one of the indicators was available. 
Because we had no information on, for example, car own- 
ership or crime rates at each area level, we could not build 
an index identical to the Townsend deprivation index or 
the Scottish Index of Multiple Deprivation used in prior 
studies. 35 However, we defined household crowding (resi- 
dential area (m 2 ) per person) for each spatial unit, as 
household overcrowding is a variable also used in the 
Townsend deprivation index. 35 A small portion of grid 
database information was missing as Statistic Finland does 
not release information on areas with <10 residents, those 
with missing data were excluded from the analyses. 

Covariates 

The covariates obtained from the employers' administra- 
tive records were age, sex and occupational title. 
Occupational tides were used as one indicator of individ- 
ual level socioeconomic status (SES). We classified indivi- 
duals into three groups: high = upper grade non-manual 
workers (eg, physicians and teachers) , intermediate=lower 
grade non-manual workers (eg, registered nurses, techni- 
cians) and low=manual workers (eg, cleaners, mainten- 
ance workers) based on the classification of occupations. 36 
This classification is determined by the activities per- 
formed in each job and education, and we have previously 
used this in our studies. 37 38 However, it may not 
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correspond to classifications used in other countries, 
because, for example, teachers in Finland are required to 
have a university degree. Other indicators for individual 
SES were level of education (high=university degree, inter- 
mediate=high school or vocational school, low=compre- 
hensive school) obtained from Statistics Finland, and 
housing tenure (owner vs other) from the Population 
Register Center. Because the study participants were 
spread all over the country (though mainly in the 
Southern and Western parts, see online supplementary 
figure 1), and because of regional variation in mortality 
rates in Finland 39 we included a four-category area-level 
variable for county in the analyses (South, West, East and 
North obtained from Statistics Finland) . 

Statistical analyses 

The associations between the two area characteristics 
and total mortality for each spatial unit were assessed 
using Cox proportional hazards regression (PHREG pro- 
cedure of SAS V.9.2) with a robust variance estimator 
that accounts for the correlation of individuals residing 
within the same spatial units. 40 We also examined the 
spatial variance using the survival package in R software 
(using frailty with gamma distribution) , which showed 
that the spatial correlation in mortality in these data is 
very small (p for frailty=0.41, variance of random effect 
0.0017). Non-significant interaction terms between loga- 
rithmically transformed follow-up time and the area 
characteristics suggested that the proportionality assump- 
tion was not violated (all p>0.05). The results are pre- 
sented as HR with 95% CI by quin tiles of area 
characteristics, with the most favourable quintile as the 
reference group. Analyses were adjusted for age, sex, 
individual-level SES variables and county. 

In further analyses we stratified the models by individual 
SES, and ran them including all spatial units of a given 
area characteristic simultaneously. As another sensitivity 
analysis, we excluded mortality that occurred within the 
first 2 years of follow-up from the data. To study the nature 
and strength of the spatial patterning of the area 



characteristics, we calculated a pooled spatial autocorrel- 
ation index over five major town areas as suggested by 
Moran (\ARIOGRAM procedure). Pooled Moran's 
indexes were used because cohort participants were scat- 
tered over the whole of Finland, but the majority (87%) 
resided within these five town areas (see online supple- 
mentary efigure 1) . The distance between the areas within 
each town was calculated as the distance between the 
bottom-left coordinates of the 250x250 m map squares. To 
visualise the spatial distribution of the area characteristics 
we generated maps using the 250x250 m map grid over a 
sample town. 

RESULTS 

During 1.45 million person-years of follow-up 3832 parti- 
cipants died. The median follow-up period was 12 years 
(interquartile range 10.0-12.0), and median time of resi- 
dence 6.4 years (see online supplementary etable 1). 
The means and SD of the area characteristics, popula- 
tion size and numbers of participants by spatial units are 
shown in table 1. Correlations between the spatial units 
for each area characteristic are provided in the web 
appendix (see online supplementary etable 2). 

For the 250x250 m area, mortality increased linearly 
(p values for trend <0.001) with increasing quintiles of 
both area characteristics (figure 1). The hazards ratio 
for mortality in the top versus bottom quintile of socio- 
economic deprivation was 1.36 (95% CI 1.21 to 1.52), 
and that of crowding 1.21 (95% CI 1.11 to 1.38) in the 
adjusted models. 

In the analyses for the lxl km map grid squares, 
crowding was linearly (p value for trend <0.001) asso- 
ciated with increases in mortality, the HR for the top 
versus bottom quintile being 1.20 (95% CI 1.05 to 1.37). 
For deprivation the association was weaker (HR 1.10, 
95% CI 1.01 to 1.20, for the top vs bottom quintile). 

When using the 10x10 km grids, zip-code areas and 
towns, household crowding was again linearly associated 
with mortality (p values for trend <0.001 10x10 km, 0.02 



Means and SD of area characteristics by spatial units and numbers of participants included in the analyses 





Spatial unit* 














250x250 m 










Variable 


N 


Mean (SD) 


1x1 km 


10x10 km 


Zip-code area 


Town 


Number of units (n) 


t 


27 495 


10 272 


1296 


1858 


422 


Unit size (km 2 ) 


t 


0.0625 (0) 


1 (0) 


100 (0) 


34.3 (96.5) 


377 (370) 


Participants per unit* (n) 


146600 


29 (32) 


227 (234) 


5777 (4260) 


762 (642) 


13073 (8883) 


Socioeconomic deprivation (z-score) 


139600 


-0.28 (0.68) 


-0.33 (0.65) 


-0.43 (0.64) 


-0.44 (0.66) 


-0.51 (0.60) 


Annual income^ (k€) 


130278 


18.5 (4.8) 


17.9 (3.7) 


17.7 (2.7) 


17.7 (3.1) 


17.4 (2.5) 


Low education (%) 


130376 


30.7 (10.9) 


32.0 (8.5) 


32.8 (5.6) 


31.8 (6.8) 


32.2 (5.0) 


Unemployment rate (%) 


139595 


8.5 (6.7) 


8.4 (4.5) 


10.8 (3.5) 


7.8 (3.4) 


7.9 (2.5) 


Household crowding (m 2 per person) 


134075 


37.7 (7.8) 


37.7 (5.0) 


37.6 (2.4) 


35.1 (2.8) 


35.3 (1.7) 



*ln statistic Finland's linkage dataset, the SD of the total number of residents in the spatial units with study participants is 285 (294) per 
250x250 m unit, 2365 (2351 ) per 1 x1 km unit, 68 251 (53 1 53) per 1 0x1 0 km unit, 8096 (51 33) per zip-code area and 1 52 355 (1 22 064) per 
town. 

fThe number of participants per unit ranges between 1 and 21 921 , depending on the spatial unit. 
^Median income per resident aged >18 with taxable income in the spatial unit. 
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Figure 1 Mortality by (A) area deprivation and (B) household crowding. HRs (95% CIs) using five alternative spatial units. 
Models are adjusted for age, sex, occupational status, level of education, housing tenure and county. 



zip code and <0.001 town). The association was the stron- 
gest in the 10x10 km grid (1.28, 95% CI 1.10 to 1.50, for 
the top vs bottom quintile) , and the magnitudes of HRs 
for zip-code areas and towns were only slighdy lower than 
those for the smaller grid based units (figure 1). 
Deprivation, however, was not associated with mortality in 
the larger units. To investigate potential confounding by 
individual-level SES, we ran the analyses for area depriv- 
ation stratified by individual-level SES variables in the smal- 
lest spatial unit. Mortality increased linearly by area 
socioeconomic deprivation in the 250x250 m square 
within each occupational group, each level of education 
and by housing tenure (figure 2). We then examined 
whether any given spatial unit drove the associations. In 
the analyses including all spatial units simultaneously we 
found the strongest associations for deprivation in the 
250x250 m area (table 2) . Associations for crowding 
became non-significant, possibly because of the high cor- 
relation between crowding in the 10 km grids and 
crowding in towns (Pearson r=0.79). Owing to this multi- 
collinearity we analysed crowding in the 250 m, 1 km and 
10 km grids and in zip-code areas simultaneously, which 
resulted in significant associations for all units except 
zip-code area (HR in the top vs bottom quintile 1.12, 95% 



CI 1.00 to 1.26 for 250 m grid, and 1.20, 95% CI 1.02 to 
1.42 for 10 km grid) . Analyses of the 250 m and 1 km grids 
with zip-code areas and towns also produced significant 
associations (HR 1.12, 95% CI 0.99 to 1.26 for 250 m grid, 
and 1.12, 95% CI 0.97 to 1.29 for town). 

When the first 2 years of mortality follow-up were 
excluded, the effect estimates for the top versus bottom 
quintile attenuated slighdy in the 250x250 m unit (1.31, 
95% CI 1.16 to 1.47 for deprivation, and 1.18, 95% CI 
1.08 to 1.28 for crowding), and in the lxl km unit (1.06, 
95% CI 0.97 to 1.15 for deprivation and 1.18, 95% CI 
1.05 to 1.33 for crowding). In the larger units, associa- 
tions between crowding and mortality remained similar 
to those in the whole data, and deprivation was not asso- 
ciated with mortality (data not shown) . 

The spatial patterning of the area characteristics, 
based on pooled Moran's Indexes for the 250x250 m 
spatial units in the five major town areas were modest — 
0.043 (95% CI 0.013 to 0.074) for socioeconomic depriv- 
ation and 0.030 (95% CI 0.019 to 0.042) for household 
crowding. Figure 3 illustrates how broadening spatial 
units affect the identification of local areas with disad- 
vantages in an example town. A substantial proportion 
of the 250x250 m squares that belong to the least 
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favourable quintile in terms of socioeconomic depriv- 
ation and household crowding may remain unidentified 
if disadvantage is defined on the basis of the surround- 
ing areas instead of the local area itself. 



DISCUSSION 

This study has two key findings. First, for socioeconomic 
deprivation, the smallest spatial unit, the 250x250 m 
square, captured the mortality associations best. The asso- 
ciations were substantially weaker when deprivation was 
defined for the lxl km squares, and no association with 
mortality was found for wider spatial units. Second, we 
found a graded association for household crowding and 
increased mortality risk across all five alternative spatial 



units. There was evidence that both small and broad 
spatial units were representative area definitions, suggest- 
ing that high household crowding in proximal areas and 
as a town average captured a partially non-overlapping set 
of mortality risk factors. Discordant findings across differ- 
ent socioeconomic exposures empirically illustrate the 
Modifiable Area Unit Problem and suggest that differ- 
ences in the spatial units used in the analyses are a source 
of heterogeneity in observed associations between resi- 
dential area characteristics and the risk of death. 

Comparison with other research 

Our findings are in agreement with several previous 
studies. An ecological study in France, for example, 
observed linear mortality effects in the two smallest area 
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Figure 2 HRs (95% CIs) for mortality by quintiles of socioeconomic deprivation in 250x250 m grids by (A) occupational status, 
(B) level of education and (C) housing tenure. Models are adjusted for age, sex and county. 
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Table 2 HRs (95% CIs) for mortality within each area 
measure in association with socioeconomic deprivation 
and household crowding when simultaneously adjusting for 
all area definitions 
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All models include age, sex, occupational status, level of education, 
housing tenure, county and simultaneous adjustments for 
'Socioeconomic deprivation for 250x250 m, 1x1 km, 10x10 km, 
zip-code area and town. 

fHousehold crowding for 250x250 m, 1x1 km, 10x10 km, zip-code 
area and town. 

scales ('commun' and 'canton') but not for the three 
larger area units. 13 Similarly, a recent meta-analysis found 
that area SES-mortality associations were slighdy stronger 
in the small (relative risk 1.10, 95% CI 1.06 to 1.15) than 
the large (1.05, 95% CI 1.04 to 1.06) area units. 18 
However, only age and sex were controlled for in these 
models. A study from Massachusetts, USA, reported 
similar effects for large and small area units — the mortality 
incidence rate ratio in the area of lowest versus highest 
socioeconomic position was 1.39 (95% CI 1.33 to 1.46) 
using the zip-code area, and 1.31 (95% CI 1.25 to 1.38) 
using the block group. 4 For crowding a significant associ- 
ation was seen at the block level (incidence rate ratio 1.43, 
95% CI 1.23 to 1.67) but not at the level of zip-code areas 
(1.18, 95% CI 0.69 to 2.00). 4 These effect estimates are in 
agreement with ours (HR 1.20, 95% CI 1.11 to 1.31 in the 



250x250 m grids, and 1.14, 95% CI 1.03 to 1.25 in the 
zip-code area) , although in our analysis the association at 
zip-code level also reached statistical significance. Stronger 
associations for the smaller of two spatial units have also 
been reported in relation to poor self-rated health, a pre- 
dictor of overall mortality. 28 

Our findings are not in agreement with data on house- 
hold income from 14 US states in which similar effect esti- 
mates for two spatial units were observed (HR for 
mortality per US$10 000 lower median tract-based house- 
hold income 1.15, 95% CI 1.13 to 1.16, and zip-code-based 
income 1.16, 95% CI 1.14 to 1.17). 11 Obviously, this com- 
parison is cruder than ours which was based on income 
quintiles. Another recent study reported linear mortality 
associations at tract-level with age-adjusted HRs in the most 
versus least deprived area — 1.53 for women and 1.66 for 
men, but when adjusted for individual-level risk factors 
these HRs attenuated considerably — 1.13 and 1.17, 
respectively. 15 In a Finnish study, mortality rate ratio for 
the highest versus the lowest proportion of manual 
workers were similar at subdistrict level (1.13, 95% CI 1.01 
to 1.25), and whole district level (1.10, 95% CI 0.98 to 
1.23). 16 At least two further studies, examining other 
health outcomes, have reported similar neighbourhood 
effects regardless of the spatial unit used. 27 29 

Why might spatial unit matter 

There are several reasons for the different area effects at 
different spatial levels. First, the scattered spatial pattern- 
ing of deprivation in the 250x250 m units in our study may 
have resulted in stronger mortality associations in the 
small units than the large spatial units. Averaging over 
larger areas that are homogenous may result in a measure 
that does not capture the local conditions relevant to 
health — the association is therefore expected to be much 
weaker for the larger areas. 33 For deprivation, local charac- 
teristics may be particularly important because they either 
exert a causal effect (eg, via psychosocial pathways and 
stress) , or because they are proxies for the individual SES, 
which is what matters causally. Nonetheless, we found no 
support for the proxy explanation because the associations 
for socioeconomic deprivation in the 250x250 m square 
remained after adjustment for individual socioeconomic 
position variables. These associations were also observed 
within each level of individual SES. 

Second, the independent mortality effects of house- 
hold crowding in local and broad spatial units suggest 
that there are likely to be multiple health-related 
mechanisms operating at different spatial levels. 
Household crowding may increase mortality risk because 
of its relation to increased transmission of diseases, such 
as respiratory and infectious diseases. 41 42 Further plaus- 
ible explanations involve neighbourhood social ties. The 
habit of smoking, for example, was found to spread 
through social ties assessed using social network ana- 
lysis. 43 This effect may be better captured within small 
than large spatial units. Household crowding at the town 
level, in turn, was only weakly related to crowding in 
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Figure 3 Illustration of the influence of area choice on detecting the local variance in area characteristics. In the maps of a 
sample town on the left, the 250x250 m squares that belong to the least favourable quintile of (A) area deprivation, and 
(B) household crowding are black, and on the right, only those squares whose neighbouring squares also belong to the worst 
quintile are black (a and b). 



smaller areas and may therefore mark increased mortal- 
ity risk through other mechanisms such as the quality of 
health-related resources, because the recruitment of 
motivated and well-trained personnel in health services 
may be easier in towns in which socioeconomic disad- 
vantage is low. 30 

Limitations 

We used GPS-coordinates to link mortality and grid-based 
area data. The advantage of using map grids is that they 
can be used for the creation and maintenance of 
population-level spatial databases such as Statistics 
Finland's Grid Database. However, the lack of coordinate- 
based sociodemographic data in many countries is one 
reason why grid-based data have rarely been used for 
defining areas in socioepidemiological studies. In the near 
future, wider use of grid databases may become possible, 
as the European Forum for Geostatistics is currently devel- 
oping guidelines for datasets and methods to link 
Population and Housing Census results from 2010 to 2011 
to a common harmonised grid. 44 It has already combined 
European population grid datasets for the reference year 
2006. 45 In this study, the used grid squares were fixed, 
whereas the administrative areas varied more in geo- 
graphic and population sizes. However, as the socio- 
economic exposure variables varied less within the large 



(eg, mean annual income in towns 17 400€, SD=2500) 
than the small units (in 250 m grids 18 500€, SD=4800), 
and because categorical exposure variables were used, we 
believe the fixed versus non-fixed area choice did not have 
a substantial influence on the results. 

At least three further methodological issues are note- 
worthy. First, although coordinate-based data provide 
many advantages for small-area research, errors in convert- 
ing addresses into coordinates may occur. 46 In this study, 
all address-to-coordinate conversions were made by the 
Population Register Center. It has reported that 90% of 
the residential building locations in Finland are correct to 
within 20 m accuracy, and that the coverage is the best in 
the city plan areas (where most participants resided) f 

Second, it has been suggested that to obtain a compre- 
hensive definition for a neighbourhood, both objective 
and subjective components of area characteristics should 
be incorporated. 48 However, we did not have subjective 
assessments of the areas, such as social ties between 
neighbours. Furthermore, we neither had data on the 
known mortality risk factors such as smoking (individual- 
level confounder), or air pollution (area-level confoun- 
der), which may vary according to the SES of the indi- 
vidual or the area, respectively. 

Third, the study population was female dominated and 
consisted of Finnish (ie, mainly Caucasian ethnicity) 
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employees from the public sector. Thus, further research 
is needed to examine whether the specific spatial patterns 
observed in this study can be generalised to populations of 
other countries, or people with different ethnic back- 
grounds or of different age and gender distribution. Our 
results might not be generalisable to the unemployed or 
to populations in which socioeconomic inequalities are 
more pronounced than those in Finland. 



CONCLUSION 

Our study has demonstrated that area effects on mortality 
may vary in different spatial units depending on the 
exposure in question. This evidence suggests that the 
choice of spatial unit for analysis is a source of heterogen- 
eity in observed associations and therefore an important 
factor in understanding area effects, interpreting previ- 
ous findings and conducting future studies. 
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